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It was recently shown ( :(uant-ph/9909074 ) that parasitic random interactions between the qubits 
in a quantum computer can induce quantum chaos and put into question the operability of a 
quantum computer. In this work I investigate whether already the interactions between the qubits 
introduced with the intention to operate the quantum computer may lead to quantum chaos. The 
analysis focuses on two well-known quantum algorithms, namely Grover's search algorithm and 
the quantum Fourier transform. I show that in both cases the same very unusual combination of 
signatures from chaotic and from integrable dynamics arises. 
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The problem of quantum chaos in quantum computers 
(QC) Q has recently attracted considerable attention af- 
ter a pioneering work by Georgeot and Shepelyansky . 
These authors pointed out that residual, uncontrolled in- 
teractions between qubits might induce quantum chaos in 
the QC if the strength of the interaction exceeds a certain 
critical level, and they argued that this might destroy 
the operability of the QC. While the model considered 
by Georgeot and Shepelyansky ||] does not describe a 
particular physical realization of a QC, and in particular 
does not allow for time dependent operating of quantum 
gates, it is sufficiently generic to mimic a quantum regis- 
ter, i.e. a (static) memory of the QC, in which a state can 
be stored and from which it should be retrievable again 
at a later time. 

It is clear that residual interactions between the qubits 
will in general lead to eigenstates of the quantum regis- 
ter that are not the original product states |0) . . . |0)|0), 
|0)...|0)|1), ... (called multi-qubit states 

in the following). Therefore, information stored in the 
register will evolve with time. If the exact eigenstates 
of the register are superpositions of only a few multi- 
qubit states, the register will oscillate quasiperiodically 
between these. However, in the case of quantum chaos, 
the eigenstates of the register are superpositions of prac- 
tically all 2" multi-qubit states. The quasipcriodic oscil- 
lations due to the interaction between qubits degradate 
then for all practical times to a decay of the original re- 
gister state into all multi-qubit states. The time scale 
for this decay is set by the inverse width of the dis- 
tribution of eigenenergies. Georgeot and Shepelyansky 
concluded that all calculations of the quantum computer 
should have finished long before the time r^^ and that 
this would limit the operability of the computer in a very 
similar fashion as decoherence. The critical interaction 
strength Jc between qubits at which quantum chaos sets 
in decreases like Jc oc 1/n with the number n of qubits, 
and the transition to quantum chaos becomes more and 
more abrupt with increasing n 

Later, Silvestrov et al. questioned Georgeot's and She- 
pelyansky's conclusions 0. They showed that even in 



the case of quantum chaos error correction schemes [HH] 
are capable of dealing with errors generated by quantum 
chaos. However, much more error correction is needed 
than in the absence of quantum chaos. Their model again 
included only static interactions between qubits. 

In this work I examine a question that is somewhat 
complementary to the one investigated in |2[|^ : Is it pos- 
sible that already the interactions introduced on purpose 
between qubits in order to operate the quantum com- 
puter lead to quantum chaos, even if there are no residual 
parasitic interactions between qubits? This is an intrigu- 
ing question, since it has implications for the amount of 
resources necessary for implementing a given algorithm. 
Already classically chaotic algorithms (e.g. for calcu- 
lating the time evolution of a classical chaotic system) 
require much more computing resources than integrable 
ones: Since errors in the initial conditions amplify ex- 
ponentially in time, one needs to spend many more bits 
than for calculating an integrable time evolution over the 
same time. 

In quantum mechanics chaos manifests itself by a sensi- 
tivity not to the initial state but to the control parameters 
^ (amongst other signatures, see below). The fidelity 
\{ip{k)\i^{k))\^ of a wave-function \tlj{k)) = U'^WW 
with respect to a wave-function \ip{k)) = U''{X)\4') de- 
creases exponentially with the discrete time k, if A and A 
are two slightly different system control parameters and 
the time evolution U{X) is chaotic. In the case of in- 
tegrable quantum dynamics the fidelity typically shows 
quasiperiodic oscillations in k. Therefore a chaotic quan- 
tum algorithm will need more resources in the form of 
more precise quantum gates, or as pointed out in 
more error correction. 

A quantum algorithm (QA) is uniquely defined by the 
unitary transformation U it induces in the entire multi- 
qubit Hilbcrt space, and the question of quantum chaos 
can be studied directly on the level of that unitary trans- 
formation, without the need to deal with the time depen- 
dent Hamiltonian by which it is generated. This is what 
I am going to do in this work, therefore being able to go 
beyond the models with time-independent Hamiltonians 
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studied in If]. 

Obviously, the answer to the question posed depends 
on the quantum algorithm. A quantum algorithm simu- 
lating a quantum chaotic system ||^,^ is by definition a 
unitary transformation showing quantum chaos, and will 
thus need very precise tuning of the control parameters. 
In this paper I focus, however, on two of the most well- 
known quantum algorithms, where the answer is less clear 
from the beginning, namely Grover's search algorithm ||] 
and the quantum Fourier transform (QFT). The latter is 
the center piece of several important QAs, like phase esti- 
mation, order-finding, the hidden subgroup problem (see 
|10|), and, most prominently, Shor's factoring algorithm 
1 11| . Also from the pure quantum chaos point of view the 
question of quantum chaos in these algorithms is very in- 
teresting. In fact it turns out that both algorithms have 
symmetry properties that lead to a remarkable and very 
non generic mixture of signatures of quantum chaos and 
quantum integrability. 

The first thing that comes into mind for checking for 
quantum chaos, is the eigenvalue and eigenvector statis- 
tics of U: It is believed that an eigenvalue- and eigen- 
vector statistics of U corresponding to Dyson's circu- 
lar ensembles indicates quantum chaos [ p"2|jr^ ]. That is, 
the eigenvalues = exp(— itpi) should show universal 
level repulsion ||l^. In the limit of large N one ex- 
pects a distribution P{s) of nearest neighbor spacings 
Si = N{ipi^i — ipi)/2Tr that is well described by the uni- 
versal Wigner-Dyson statistics [|l5| . If C/ is covariant un- 
der any anti-unitary operation T that squares to unity 
there is always a basis in which the eigenvectors of U can 
be chosen real. The relevant random matrix ensemble 
is then the circular orthogonal ensemble (COE) with a 
P{s) very well approximated by 



the searched one in which case the oracle changes the 
state's phase by tt. Grover's algorithm commences with 
the Hadamard transformation 



where 
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is a Hadamard transformation for the ith qubit. Starting 
from the register state |0)...|0), the Hadamard trans- 
formation brings the system into a superposition of all 
multi-qubit states with equal weight 1 / VN. Next comes 
an iteration of the oracle query followed by a "diffu- 
sion matrix" D. The latter has matrix elements Dij = 
2/N — Sij. The total algorithm thus reads 



Ug = [DOyH . 



(5) 



The optimal value for the integer p is given hy p — 
|7r/(40)] with sin^ 9 — 1/N ([.] denotes the integer value) 
||l7| , and I have chosen this value for all calculations. Re- 
markably, the quantum computer finds the searched ele- 
ment with ^ \/N queries, whereas a classical computer 
when presented with the same problem would have to 
ask the oracle on the average N/2 times. 

From the definition of Uq it is clear that C/g is real. 
Covariance under conventional time reversal symmetry. 
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thus is fulfilled if Uq is symmetric, Ug = Uq, where Uq 
denotes the transposed matrix. Using the commutation 
relations between H, O and D (or by evaluating the ma- 
trix numerically, see fig.|l|) one convinces oneself that Ug 
is almost symmetric. 



The best known example is conventional time reversal 
symmetry, in which case T is the complex conjugation 
operator @. The eigenvector statistics is usually de- 
scribed in terms of a distribution of eigenvector com- 
ponents. Picking any component Ci of any eigenvector 
at random, random matrix theory (RMT) predicts for 
y — N\ci\'^ in the COE case the so-called Porter-Thomas 
distribution, 
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Let us now have a look at Grover's algorithm It 
allows to find an entry with index ^ in a unsortcd quan- 
tum database that is distinguished from the others by a 
given property. The distinction may be formalized by an 
oracle query O which in the multi-qubit basis is a uni- 
tary diagonal matrix with entries On = 1 — 25i^ where Si^ 
is the Kronecker delta. Thus, presented a register state 
the oracle always gives back the same state unless it is 
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FIG. 1. Average absolute value of the matrix elements of 
the symmetric part of Uq (circles) and the antisymmetric 
part (squares). The former decays like 1/y/N, the latter like 
where A'^ — dimn is the dimension of the Hilbert space. 
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FIG. 2. Eigenphase spectrum for the Grover algorithm for 
q = 7 qubits, ^ = 64. The spectrum is symmetric and — tt 
and n should be identified. Only two independent eigenvalues 
differ from ±1. 

In fact the average absolute value of the matrix ele- 
ments of the symmetric part of Ug decays like 
whereas the average absolute value of the antisymmet- 
ric part decays like Just one element escapes from 
the general rule set by the average: Ug^o, the element 
in the zeroth column pertaining to the searched element 
is of order unity minus a 1/\/N correction (that is how 
the whole algorithm works), whereas Ugq^ plays no such 
special role. Nevertheless, for ^ oo the weight of 
this particular element is zero and therefore Ug has the 
unusual property that it is at the same time unitary as 
well as real and (almost) symmetric. Thus, all eigenval- 
ues have to be at the same time on the unit circle, and 
have to be real. Therefore all eigenvalues are expected 
to be either unity or minus unity! This is well confirmed 
numerically (see FIGJ|). In fact it turns out that all 
eigenvalues, even those which are not plus or minus one, 
are to very good approximation 6th roots of one, and 
Grover's algorithm is therefore approximately a 6th root 
of the unit operation! 

The resulting high degeneracy of the eigenphases is 
in strong contrast to the level repulsion that goes along 
with quantum chaos. However, absence of level repulsion 
in U does not exclude quantum chaos. Indeed, suppose 
some U did show level repulsion, then C/*^ with M ^ I 
will in general not, since the spectrum gets completely 
mixed by winding it AI times around the unit circle. A 
Poissonian statistics will be the consequence. Neverthe- 
less, if U has a chaotic classical counterpart, so will U^'^ , 
and to conclude that U^^ is not chaotic just from the 
level statistics is therefore not possible. If the phases of 
U are commensurate one can even find an M so that all 
eigenvalues are degenerate. 
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FIG. 3. Integrated eigenvector distribution 
I{y) = R{x) dx for the Grover algorithm for n = 7 and 
various values of ^. The different symbols denote different 
selected numbers (^ — 0, 3, 5, 64 and 127). The full line is the 
integrated Porter Thomas distribution. 

Let us look at the eigenvector statistics of Ug- FIG.^ 
shows that is obeys almost perfectly RMT! This is, how- 
ever, a pure consequence of the high degeneracy. Indeed, 
any linear combination of eigenstates in the subspace per- 
taining to the degenerate eigenvalue A = 1 is again an 
eigenstate to the same eigenvalue (and correspondingly 
for A = —1). Thus, a set of eigenvectors can be oriented 
in a completely arbitrary way in this subspace with the 
only restriction that they be all normalized and mutually 
orthonormal. This is just as in RMT, where the eigenvec- 
tors arc all statistically independent from the eigenvalues, 
and their orientation is statistically uniform on the whole 
possible hyper-sphere. In the present situation the diag- 
onalization routine picks initial orientations at random 
and thus mimics perfect RMT behavior in both of the 
two subspaces. 

Given the unclear picture presented by the eigenvalue 
and eigenvector statistics, it seems reasonable to examine 
directly the sensitivity with respect to slight variations in 
Ug- One might do so by applying Peres' original scheme, 
i.e. creating one perturbed algorithm Uq and studying 
the decay of overlap between Uq\iP) and (C/q)"|'0) as 
function of n for some random initial {ip}- A natural 
parameter in Ug that might be varied is the number of 
iterations p of the transformation DO. Forp~ \/]V going 
from p to p + 1 appears indeed as a small perturbation. 
However, all that has been said about the spectrum of Ug 
applies to the such perturbed Uq as well, i.e. all eigenval- 
ues will generically be ±1. A spectral decomposition of 

Ug = \u^'>){u^^\-J:i l"^^>(«l"VE.e''^H^fc)KI 
(where (pk are the phases different from and ±7r) and 

correspondingly for Uq shows that for even n the over- 
lap {^\{Uq)"{Uq)'^\'4j) depends on n only due to the ex- 



ceptional eigenvalues that differ from ±1. All the other 
eigenvalues contribute only two different values to the 
overlaps, depending on whether n is even or not. There- 
fore, varying p (or introducing any other perturbation 
that does not lift the degeneracy of the spectrum) does 
not lead to exponentially decaying overlap. 

I have therefore applied another perturbation: All DO 
factors in Ug are multiplied with a random matrix Vi 
(drawn independently for each factor) close to unity, 

= iDO)Vi . . . {DO)VpH . (7) 

I constructed all of these random matrices Vi as ten- 
sor products of 2 X 2 orthogonal matrices 0{2,ip) close 
to diag(l, 1) in the Hilbert space of each qubit, Vi = 
0(2, (fii) (g) . . . (g) 0(2, fpi), where fu, . . . , fpi are chosen 
randomly and independently from a uniform distribution 
in the interval — e/2 . . . e/2, and 0(2, </?) is a 2 x 2 orthog- 
onal transformation acting on a single qubit as 

0{2,^)^( ^^"^V (8) 

\^ — sm If cos ip J 

Fig.^ shows quasiperiodic overlap with almost perfect re- 
vival within 500 iterations for a 5 qubit Grover algorithm 
with ^ = 2 as selected index and e = 0.1. For iteration 
numbers exceeding 5000 one notices a slight decay, but 
nevertheless the Fourier spectrum of the decay shows a 
few very sharp and strong peaks, indicating the quasi- 
periodic nature of the function. We conclude that ac- 
cording to Peres' criterion, Grover's search algorithm is 
free of quantum chaos. 

One might object that as a quantum algorithm, Ug will 
typically not be iterated, though. I therefore also applied 
another criterion proposed by Schack et al. which 
seems to be more appropriate in the current situation. 
These authors examined the distribution of angles be- 
tween vectors propagated by slightly perturbed unitary 
matrices, and embedded Peres' original sensitivity cri- 
terion into an information theoretical framework. They 
showed for the example of a kicked top that in the case 
of a chaotic quantum map the distribution P{a) of an- 
gles a between Hilbert space vectors propagated by many 
slightly and randomly perturbed unitary transformations 
corresponds to that of randomly chosen vectors, which re- 
sembles a Gaussian peak. An overall average angle can be 
steered with a deterministic part of the random vectors 
and adapted to P{a) of the Hilbert space vectors. On 
the other hand, for integrable quantum maps the random 
perturbations lead to many more or less degenerate an- 
gles, as the propagated vectors do not explore all Hilbert 
space dimensions. The angle distribution therefore typ- 
ically contains several more or less pronounced peaks. 
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FIG. 4. The overlap between a state propagated by 
Grover's algorithm and by a slightly perturbed Grover al- 
gorithm (q = 5 bits, (, — 2, e = 0.1) plotted on a logarithmic 
scale. The overlap shows quasiperiodic oscillations, charac- 
teristic of integrable behavior. 

Instead of adapting the deterministic part of the ran- 
dom vectors, I chose to "unfold" the angles, i.e. to rescale 
them by the average angle. It turned out that "univer- 
sal" distributions are obtained by this procedure, both 
for the randomly drawn vectors as well as those propa- 
gated by Uq. In the latter case one obtains a distribution 
which over several orders of magnitude of the perturba- 
tion strength e is independent of e, ^ in the former case 
the unfolded angle distribution is independent of the de- 
terministic part of the random vectors. 

I have constructed two particular classes of pertur- 
bations. The first one closely corresponds to the orig- 
inal recipe by Schack et al. ||l^. The perturbed algo- 
rithm is obtained by multiplying all DO factors in Ug 
with only one out of two random orthogonal matrices 
close to unity, namely V+ = 0{2,ipi) (g) . . . ,0{2,Lpp) 
or V- = V^^, where again ipi,...,ipp are chosen 
randomly and independently from a uniform distribu- 
tion in the interval —e/2... e/2, but are kept fixed for 
all factors; i.e. we obtain 2^ perturbed Grover algo- 
rithms DOV+DOV+ . . . DOV+, DOV+DOV+ . . . DOV^, 

. . ., DOV-DOV DOV-. A random initial vector is 

then propagated by these 2^ perturbed matrices and the 
angle distribution between the resulting vectors is ana- 
lyzed. 



^This is not true, though, for the perturbed Grover algo- 
rithm using the 'digital' form of perturbation; see below. 
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Fig. H shows that this sort of 'digital' perturba- 
tion does not entirely randomize the propagated vec- 
tor. Rather the distribution of angles shows many pro- 
nounced peaks which means that different perturbations 
lead to similar Hilbert space vectors. Thus, from this 
plot one would conclude that the behavior is integrable. 
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FIG. 5. The distribution of angles between Hilbert space 
vectors propagated with 2*" perturbed Grover algorithms 
(g = 5, ^ = 2, e = 0.1 full line) is compared to the distri- 
bution of angles between random vectors (dashed line). In 
both cases the angles were "unfolded" such that the average 
angles equal unity. 
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FIG. 6. The distribution of angles between Hilbert space 
vectors propagated with 50 perturbed Grover algorithms 
{q = 5, — 2, e — 0.01 ; independent random matrices 
for each factor DO, full line) is compared to the distribu- 
tion of angles between random vectors (narrow peak, dotted 
line). All angles were "unfolded" such that the average angle 
is unity. Also shown is the corresponding distribution for the 
QFT (g = 5, e = 0.01, dashed line). 



However, the situation is very different for the second 
form of perturbation which corresponds to the one intro- 
duced for studying the decaying overlap, i.e. instead of 
choosing between just two matrices V- and V+ each fac- 
tor DO is multiplied with an independent random matrix 
of the form V{ipu, . . . , ippt) = 0(2, ipu) (g) . . . (g) 0(2, ippi). 
Fig. 1^ shows that in this case a broad, Gaussian like dis- 
tribution without much further structure arises, much 
broader in fact (after unfolding) than the angle distribu- 
tion obtained from the random vectors. 

We may therefore conclude that for certain pertur- 
bations and according to the criterion by Schack et 
al. Grover's search algorithm does show hypersensitivity 
with respect to perturbations. This is quite surprising in 
the light of Grover's earlier finding that almost any uni- 
tary transformation substituted for the Hadamard ma- 
trix still leads to a functioning search algorithm , and 
it seems, indeed, that the criterion by Schack et al. is sin- 
gled out compared to the other criteria examined. Note, 
however, that there is not necessarily a contradiction to 
Grover's finding, since for Grover's algorithm to work it 
is enough that only the first column of Uq be largely 
unaffected by the perturbations (and the same argument 
applies to all quantum algorithms that start from just one 
initial state, typically the state |00...0))! The present 
result is more general in as much as it makes a statement 
about the sensitivity of the entire matrix with respect to 
perturbations. 

How rapidly does the average fidelity decrease with 
the perturbation strength? If we measure lack of fidelity 
as average absolute error of all matrix elements, the an- 
swer is: linearly, over several orders of magnitude of e 
(0.0001 < e < 0.1), as can be seen from Fig.0, and as 
would have been expected from perturbation theory. 

Let me now come to the quantum Fourier transform 
(QFT). Since it is a universal part of several proposed 
quantum algorithms including Shor's algorithm ||ll|] , it 
makes sense to give the QFT special attention. The 
quantum Fourier transform Uft on a n bit register (with 
qubits indexed as ... n — 1) can be constructed from 
one- and two-qubit operations as ||ll| 

Uft = FHoSo^iSq^2 ■ ■ ■ <5'o,ri-l^^l • • ■ Hn-3Sn-3,n-2 

Hn-1 , (9) 

where Sj,k is a conditional phase shift matrix between 
qubits j,k defined by Sj_k — diag(l, 1, 1, exp(i7r/2'^~-' )) 
in the basis 00, 01, 10, and 11 formed by the two qubits 
j and k. The matrix F flips all qubits. With all the two- 
qubit interactions introduced, one would naively expect 
quantum chaos. However, the QFT is constructed such 
that in the whole Hilbert space C/pT is very simple and 
symmetric. 



UpTik 



1 



^i2Trlk/N 



N 



(10) 



5 



One easily convinces oneself that = l! Thus, 
all possible eigenphases are 0, tt, and ±7r/2. The 
situation is therefore even simpler than in Grover's 
algorithm: There is an exact relation that dic- 
tates a high degeneracy of only four possible eigen- 
phases. So again, there is no level repulsion. 
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FIG. 7. The absolute average error of the matrix elements 
of Grover's search algorithm Ug (in units of the absolute aver- 
age matrix element) for random perturbations (independent 
matrices V{ip) for each factor DO); q = 5, a — 2. 



The matrix C/pT is covariant under conventional 
time-reversal: KUftK^^ = U* = W since U = 
U"^ . And a numerical evaluation of the eigen- 
vector statistics leads again to good agreement 
with the Porter-Thomas distribution, which, how- 
ever, results once more from the high degeneracy 
of the eigenvalues and not from quantum chaos. 
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FIG. 8. The overlap between a state propagated by the 
QFT and by a slightly perturbed QFT {q = 5 bits, e = 0.1) 
plotted on a logarithmic scale. The overlap shows quasiperi- 
odic oscillations, characteristic of integrable behavior. 

For studying the sensitivity of U-pT with respect to per- 
turbations, I slightly perturbed all phases in the condi- 
tional phase gates by adding an additional random phase 
with the same (on the average) relative amount for all 
gates . Fig. ^ shows that the overlap of a random state 
propagated with a slightly perturbed matrix C/prp and the 
same state propagated with the original Uft shows again 
quasiperiodic oscillations. The basic period is two, as 
could be expected from (C/|rp)jj — Si^jy-j, i-c the exact 
QFT leads back to the same starting vector up to a rela- 
beling of the indices when applied twice. The perturbed 
QFT only leads to a small deviation from this relabeled 
initial vector, and therefore almost perfect overlap is re- 
stored every second iteration. Thus, according to Peres' 
criterion the QFT is free of quantum chaos. 

Fig. ^ shows the angle distribution resulting from 
propagating a random initial vector by 100 differ- 
ently perturbed matrices for an algorithm running 
on five qubits. The Hilbert space vectors are com- 
plex now, so I defined the angles ai^k as ai^k = 
arccos ( |( V'i I V'fe ) I / \/ ( V'i I V'i ) ( V'fc I V'fc ) ) , and aU angles were 
again rescaled according to ai^k — > on^kl {ct), where (a) 
is the average angle for all pairs i ^ k. The resulting 
distribution depends for a small number q of qubits still 
substantially on q. For q = 3 a distribution is obtained 
that is remarkably close to the famous Wigner Dyson 
surmise for the orthogonal ensemble (^, but for larger 
numbers of qubits the distribution approaches more and 
more a Gaussian, not too different from the distribu- 
tion obtained from the Grover algorithm (see Fig. |^). 
Again the distribution is stable in the range I examined 
it, 0.0001 < e < 0.1 . Thus, also the QFT shows hy- 
persensitivity with respect to random perturbations as 
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judged by the criterion of Schack et al., and in contrast 
to Peres' original criterion! 

On the other hand it is known that if the controUed 
phase gates are performed to a precision A = l/p{n) 
where p is a polynomial in the number of qubits, then the 
maximum error of the final state Uqft\iP) for all input 
states is of order -n? /p{n) [Q. Thus, only polynomial 
precision is needed for the controlled phase gates. Cop- 
persmith has shown, indeed, that an approximate QFT 
can be obtained by dropping the gates with the exponen- 
tially small phases altogether . 

In summary I have shown that both Grover's search 
algorithm and the QFT give rise to the same unusual 
combination of quantum signatures of chaos and of inte- 
grability. Strong symmetries lead to a large degeneracy 
of the spectra of eigenvalues of the unitary matrices rep- 
resenting these algorithms. In fact, the QFT is a fourth 
root of the unity matrix, and Grover's algorithm is to a 
good approximation a 6th root of unity! The correspond- 
ing lack of level repulsion would be commonly interpreted 
as absence of quantum chaos. The eigenvector statistics 
closely follows RMT predictions, and one would com- 
monly interprete this as a signature of quantum chaos. 
However, as shown above, it is here but an artefact aris- 
ing from the highly degenerate spectrum. The overlap 
between a random state propagated by a perturbed algo- 
rithm and the same state propagated by the correspond- 
ing unperturbed algorithm shows quasiperiodic oscilla- 
tions both for Grover's algorithm and the QFT, thus 
signaling absence of quantum chaos, in agreement with 
earlier studies addressing the stability of these codes. 
The only criterion indicating quantum chaos and not ev- 
idently explainable by an artefact is the distribution of 
angles between vectors propagated once by many slightly 
disturbed algorithms. After unfolding the angles univer- 
sal distributions are obtained for large enough Hilbert 
spaces, both for Grover's algorithm and for the QFT, 
that resemble the one for random vectors. 
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